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We summarize recent developments in the field of higher dimensional bosonization made by the 
£ — ' authors and collaborators and propose a general formula for the field operator in terms of currents 

and densities in one dimension using a new ingredient known as a 'singular complex number'. 
Using this formalism, we compute the Green function of the homogeneous electron gas in one 
spatial dimension with short-range interaction leading to the Luttinger liquid and also with long- 
( — i ' range interactions that leads to a Wigner crystal whose momentum distribution computed recently 

exhibits essential singularities. We generalize the formalism to finite temperature by combining with 
the author's hydrodynamic approach. The one-particle Green function of this system with essential 
singularities cannot be easily computed using the traditional approach to bosonization which involves 
the introduction of momentum cutoffs, hence the more general approach of the present formalism 
is proposed as a suitable alternative. 
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As the term suggests, 'bosonization' is an effort to recast theories involving entities that are not bosons in terms 
of bosons that are typically expressed in terms of bilinears of the original fields. Indeed it is even possible to recast 
theories involving bosons in terms of other bosons. This activity is not merely a pedantic exercise, for the end 
£Zj ■ result of these efforts is a nonperturbative technique for studying the original theory. In particular, applying this 
technique to fermions in one spatial dimension leads to what are known as non (Landau)-Fermi liquids where the 
i momentum distribution is continuous at the Fermi momentum. Thus one can bosonize fermions, spins and even 

complex scalar fields which are themselves bosons. Bosonization, according to our understanding, is nothing but the 
t-H , polar representation of a complex number. Bosonization of spins is accomplished by polar decomposing the ladder 
operators which leads to a semiclassical theory of spins. Complex scalar fields may be easily bosonized by polar 
decomposition as well. One may suspect that a similar decomposition should be feasible for fermions too. This 
has proved harder than one might hope. In one spatial dimension, the Thirring model which describes relativistic 
fermions self-interacting via a short-range repulsion was shown to be equivalent to the so-called Sine-Gordon theory 
which involves scalar fields 1 -. This technique has been adapted to study condensed matter problems where the parabolic 
[~^. ' dispersion of the free fermions in Galilean invariant systems is linearized near the Fermi points so that the fermions 
now have linear dispersion and are moving with the Fermi velocity. The Fermi velocity takes on the role of the speed 
of light making the analogy complete. In more than one dimension, this program has not been particularly successful. 
In hindsight, it appears that it is premature to regard the framework available even in one spatial dimension as a 
closed subject. The Luttinger model caricature it seems, is unable to handle some exotic situations involving truly 
long range interactions in condensed matter systems. This was pointed out by the present author and this work 
is a summary and continuation of earlier work o 23456 . The work of Schulz 7 on electrons interacting with potential 
(V(x) ~ using standard bosonization techniques is able to study the structure factor the Wigner crystal quite 

effectively but when it comes to the one-particle Green function, the results are very sketchy and cutoff dependent. 
It is not clear whether the method of Schulz can be used to study the one-particle Green function of a ID Fermi 
• i-H , system with potential V{x) ~ — |x|, since our results show that the momentum distribution does not depend on 
arbitrary cutoffs, whereas they are mandatory in such standard bosonization methods. Furthermore, there are some 
technical subtleties involving so-called Klein factors that ensure fermion commutation rules between Fermi fields in 
the traditional approach that are far from satisfactory. The purpose of the present article is to summarize recent 
developments in the subject made by the author and his collaborators and to present a formula for the field operator 
in terms of currents and densities that is valid in a general sense. We stress that this write up is by no means a 
review of the entire field of bosonization, it only highlights the author's contributions, while not completely ignoring 
the work of others. We then go on to use this formula to calculate the Green functions of interacting systems in one 
spatial dimension, one with short range repulsion leading to the Luttinger liquid and also with a specific long range 
interaction which leads to the Wigner crystal. 

At this stage it is appropriate to survey some relevant literature on the subject of bosonization of fermions in 
general and higher dimensional bosonization in particular. This subject, as is well known, started with the papers 
of Tomonaga and Luttinger. Later on, Lieb and Mattis, Luther, Luther and Peschel developed it further. All these 
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well-known literature may be found in the reviews and texts&£. This program was taken up by Haldanei^ who coined 
the term 'Luttinger liquid' to describe systems whose Green functions have power law singularities rather than simple 
poles. The generalization of these ideas to higher dimensions was attempted by Luther. The effort to generalize these 
ideas to higher dimensions seemed phenomenologically ill-founded as one does not expect to see non-Fermi liquids in 
higher dimensions except perhaps for very long-range interactions which are un-physical. However the phenomenon 
of high temperature superconductivity suggested the need for such a ground state. In the mid-nineties starting from 
the work of HaldaneAA, Castro-Neto and Fradkin 12 , Houghton et. alJ^, Kopietz and collaborators^ attempted to 
develop a theory in higher dimensions along the lines of the Tomonaga-Luttinger theory. The works of Kopietz and 
collaborators deserve special mention since they have persisted with this technique over the years. These efforts 
closely mimic the theory in one dimension that is valid only in the sense of the random phase approximation (RPA)i£. 
Attempts to go beyond this approximation using the original formalism seems futile. Hence Setlur and collaborators 
developed a new scheme loosely based on the work of Castro-Neto and Fradkin 12 to overcome these difficulties. In 
what follows, we describe the outcome of these efforts made by the present author and his collaborators and go on to 
write down a formula for the field operator in terms of currents and densities. 

II. REVIEW OF THE GENERAL FORMALISM 

In what follows we describe briefly the main results in the authors' earlier works. It must be stressed that this 
write up is not a substitute for a reading of those works. Let c k and be fermion annihilation and creation 
operators. We may define new operators using these that are called sea-bosons. They are in general, complicated 
non-local combination of number conserving products of Fermi fields^. However in the sense of the random phase 
approximation(RPA) we may write, 

A k (q) « n F (k - q/2)(l - n F (k + q/2)) 4_ q/2 c k+q/2 (1) 

Here uf(p) — 6(kF — |p|) is the momentum distribution of free fermions. This object A k (q) has been shown to obey 
the following commutation rulesi 

L4 k (q),4,(q')] =n F (k-q/2)(l-n i ,(k + q/2)) <W<W I [A k (q), A* (q)] = ; K(q), A{, (q )] = (2) 

The defining equation for A k (q) may be partially inverted and a formula for the number conserving product of two 
Fermi fields may be written down in terms of these sea-bosons. Define c k ,< = npik) c k and c k > = (1 — np(\i)) c k . 
At the level of RPA we may write, 4_ q/2 < c k+q / 2 ,> « A k (q) and cj c + q /2,< c k-q/2,< ~ and cj c + q /2,> c k-q/2,> ~ 0. 
Thus at the RPA level c k+q / 2 c k - q /2 ~ ^4 k (— q) + ^(q)- In general (beyond RPA) we may write^ (q 7^ 0), 

4 +q /2 c k-q/2 « ^k(-q)+4(q)+5Z4 +q /2- qi / 2 ( c li)^k-q 1 /2(-q+qi)-X) A k- q /2 +qi /2( ( li)^k +qi /2(-q+qi) (3) 

qi qi 

The above identification together with Eq.@ can be shown to be sufficient to reproduce the exact commutators 
between Fermi bilinears. The correspondence in Eq.Q is also sufficient to reproduce all the dynamical correlation 
functions of the operator cj c+q / 2 c k-q/2 °f the free theory provided we set the kinetic energy operator to be K = 

S kq ^ik( c l)^k(q) [for more details please consult our published works]. Thus any theory involving fermions 

that conserves their total number may be reexpressed in terms of these bosons. The main purpose of the present 
article is to fully invert the defining equation for A k (q) and express the field operator c p alone in terms of these bosons. 
To accomplish this, we first attempt to polar decompose the field operator in real space ip(r) — ^7= J2 P e * p r °p- 

ijj(f) = e^P^e-^sf^f) (4) 

Here we have introduced a new variable n(r) which is a canonical conjug ate to p(r) = ^ (r)i/}(r)iS.. This means that 
[n(f),II(f')] = and [n(f), p(r)] = iS(f — r) and of course [p(r),p(r)] = 0. As an operator we know that p is 
non-negative. It is then well known that n cannot be self-adjoint. In fact we may write n(r) = Xq + H(r) where H(r) 
is strictly self-adjoint and Xq is conjugate to the total number [Xo,iV] = i. In what follows we regard the number 
operator to be equal to infinity (thermodynamic limit) in which case Xq should be arbitrarily small in order for it to 
be a conjugate to the number operator. Hence we shall not be too careful and we shall treat Xq as a real c-number 
and ignore it altogether. 
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We now divert our attention and examine the property of the fermion current (density) operator 3(f) — 
Im[ip' (f)'Vip(f)]. In two and three dimensions, we may construct the operator W = p(V x J) — Vp x J. We 
first represent currents and densities in real space p(f) = ^ S(f— ft) and 3(f) — J2i S(f— r,) + J2i S(f~ n) ■ 
By acting this on fermionic wavefunctions in real space we conclude that W = 0. This means that W = p 2 V X - 3 = 0. 
In other words, there exists a scalar II such that J — —pVIL. Independently we may argue that a more general ansatz 
J = — pVII + C([p];f) fails to reproduce the current-current commutator unless C = 0. It appears then, that the 
function A in Eq.([J| should be independent of f. This A is crucial since it determines the statistics of the field i/j. 
In particular, setting A = describes bosons rather than fermions. In our earlier work£, we argued that the ansatz 
in Eq.Q may be used to derive an action in terms of II and p. In the lagrangian formalism, there are no operators. 
We demote the operators IT and p to the status of real numbers and use Eq.Q in the action for free fermions, 
S = J 1/3 dt J d d x ip^ (idt + ^)ip- This led to the following action for free fermions, 

r-w r , ( p(vn) 2 + ^ , 

S = / dt d d x\ pd t U - Vf([p];x) - '- (5) 
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Here Vf is a functional of the density that has to be fixed by making contact with the properties of the free the- 
ory. We have shown that the RPA limit of the above action may be rigorously derived using sea-bosons^. We have 
also shown£ that this leads to the following expression for the generating function of density correlations of a ho- 
mogeneous electron gas in terms of the corresponding quantity for the free theory. If Z([U}) is the function that 
generates density correlations of the homogeneous electron gas where the electrons interact with a potential t> q and 

Zo([?7]) is the corresponding quantity of the free theory ( for example (^ju^— — ^ =< Pq,nP-q,-n > and 
< Tp(q, t)p(— q, t ) >= Yln e ~ W " ' < Pq,nP-q,-n > and w n — 27m//3 is the Matsubara frequency ) then, 

Z([U}) = f D[U] e^^^-CX^-n-^-J Z ([U'}) (6) 



It appears that we have to judiciously combine the hamiltonian or the operator version and the lagrangian version in 
order to obtain useful results. Thus we wish to now revert to the operator description to try and express the field 
operator explicitly in terms of currents and densities. We observed that current-current commutator implies that A is 
independent of f. Unfortunately, this conflicts with the requirement that ip obey fermion commutation rules. Indeed, 
imposing these rules on ip in Eq.((3]) leads to the following constraint on A. 

e iA([p];r) e iA([{p(x)+8(x-r)}];r' ) _ _ e *A([p] ;r ) e iA([{p(x)+S(x-r' )}] ;F) ^ 

If A is independent of f this is impossible. Hence we seem to have reached an impasse. There is a way out this difficulty 
using the notion of what may be called 'singular complex numbers'. We describe this concept in the following section. 

III. FIELD OPERATOR USING SINGULAR COMPLEX NUMBERS 

In our earlier work£ we argued that the field operator in momentum space may be expressed directly in terms of 
the sea-bosons provided we invoke the concept of a singular complex number : w p — e~ lN where N° — > oo and £ p 
is arbitrary. Therefore, w p w p > = w p >w p — 5 p '. This rather unusual quantity may be motivated using the following 
argument. Consider the number operator nk- We may introduce formally a conjugate Pk namely an operator that 
obeys [Pk, P^ 1 ] = and [Pk, n k '] = i <5 k k ' . If we are going to treat rik as a c-number namely, = np(k) 1, then we 
have to ensure that Pk is a formally infinite c-number in order that it is a conjugate to rik- Thus u>k = e~ tPk has the 
properties that we have just described. We argue that the term e lA may be rewritten using these complex numbers. 
Let us invoke the following ansatz that is inspired from our early work 2 - : 

i)(f) = ( _L£V E »(M^ w p e lkFpp n F (p)) e- m ^ Jplf) (8) 



t/,t(f) = ^fplf) e m ^ ( -j= e-^Mtf w p e~ lkF P- p n F (p) ) (9) 
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Here E p ([p];r) is real. The idea is that the rapidly varying part is written separately as a multiplying exponent 
e ±ik F p-r^ rp^ e s ]^ ow jy var yi n g portion is in the density and phase variables. Notice that the theory presented here 
is very general and there are no momentum cutoffs at the outset. Therefore the terms 'slow' and 'fast' are merely 
suggestive of approximations that will have to be used in the practical computations where such a distinction acquires 
concrete meaning. If we postulate that VE p ([p];f) = — VE_ p ([p\; r) then we find that, p(r) = i/A (f)i/j(r) and 
J(r) = Im[ip' (r)(Vip(r))] — — p(r)Vn(r). Therefore, current algebra - the mutual commutation rules between 
currents and densities is then trivially obeyed. In other words, e iA ^ p ] ;r ^ = ^= ^2 p e l - Ep ^ p l ;r ^ w p e lkF ^' r uf{p) is 
unitary. Now we have to apply fermion commutation rules. We expect to derive a recursion relation similar to Eq.©. 
Indeed we find that both the requirement {ip(r),ip(r )} = and {tp(f), ip^{r)} = S(f— r ) are obeyed if we ensure 
that the following recursion holds, 

e -tE p ([{p(x)+8(x-f' )]; f) e iEp(b]; r) = _ e -iE p ,([{p(x)+S(r-x)}]: r ' ) jEj ([p]; r' ) ^ 

In one dimension, we try the following ansatz. Note that p = ±1 in one dimension. E p i([p]; r ) = 

J^^dx p(x) Dp/(x,r ). The recursion relation implies, e~ l Df> ^ r — — e 1 D t>'^ r ' r ■* whereas the unitarity condi- 
tion forces us to choose, Dpi (x, r ) = — p D re d(x, r ). Both these are satisfied by the choice Dpi (x, r ) = — np 6(x — r ) 
where 8(x) is the Heaviside step function ( 6(x < 0) = and 8(x > 0) = 1 ). We have to verify that this choice 
reproduces the Green function of the noninteracting theory. This is shown later on. In higher dimensions we similarly 

try, E p ([p\; f] = J d d x p(x) Dp(x,r). The recursion condition in higher dimensions : e -4 Df,( - r ,r ) — — e 1 D p'^ r ' r ^ 
seems rather hard to satisfy since p and p can be parallel to each other, antiparallel, or anything in between. If we are 
willing to require that only {ip{r), ijj'(r )} = d(r — r ) be obeyed then a simple choice might be sufficient. In any event 
it is the propagator that we are interested in and for this it is this rule that is important. In this case p — p and we 
may choose, Dp(r, r ) = tt 8(p ■ (f— r )). This choice obeys both the unitarity condition : WDp(f, r ) = -VD-p(f, r) 
and the recursion. The recursion breaks down when p ■ (r*— r ) = but these are a set of points of measure zero. 
Unfortunately, it can be shown that this choice fails to reproduce the free propagator. This suggests that the function 
E p may be a nonlinear function of the density in higher dimensions. Therefore in higher dimensions it is better to 
try and express the field variable directly in terms of sea-bosons. A partially correct formula has been provided in 
an earlier worl*£, it appears that more work is needed to make that agenda practically useful. For now we use the 
expression for the field in terms of the hydrodynamic variables to calculate the Green function in one dimension, first 
of the Luttinger liquid. In the case of the Wigner crystal we find that a reinterpretation is needed that essentially 
involves re-expressing the field directly in terms of sea-bosons. For example we may rewrite the field in one dimension 
as follows. 

^( fj t) m ( e ~ !9 "l-fJp,«-^-,(t)] Wp e ik F pr np ^ j (11) 



p 



The rapidly varying e~ lept has to put in by hand, however it may be motivated by realizing that the X q= o term is 
conjugate to the total number of particles and this picks up a contribution similar to the one suggested upon time 
evolution with respect to the hamiltonian of the system since, for example, the hamiltonian of the free Fermi theory 
may be written as H — Neo + J^k q ^A\.(q)Ak(q)- If we set eo = ep then we recover the factor suggested. We write 
ftr = Efc[4t(-9)+4fa)] =Pr(q)+Pl(q), where Pn(q) = ?2 k>0 [A k (-q)+Al(q)] and p L {q) = *£ k<0 [A k (-q) +A{{q)] 
and X q = i^(p R (-q) - p L (-q)). If p = +1 then [-^p q - iX_ q ] = -^p R {q). If p = -1 we have [£p q - iX_ q ] = 
^ZPl(q)- Thus we have the familiar result for right and left movers. We may write ip(r, t) = ip_n(r, t) + ^(r, t) or, 

^t pR{q) w p e lkFr n F (p)) + I -L^ e i:. "^^(9) w p e~ lkFT 





As usual we have [p R (q), p R (-q)] = T,k>o\. A k(-Q),A(-9)] +T,k>oi A k( ( l); A k( < l)] = ~1§ and [pdq), PL {-q)\ = §§■ 
This expression is well-known to the traditional bosonizing community and is known to reproduce the exponents of the 

Luttinger liquid correctly. In passing, we note that the present formalism does not allow terms such as tp\ { (x)^L(x ), 
this being identically zero due to the singular complex number w p . Thus in our formalism, the terms responsible 
for backward scattering come from the quadratic corrections on the right hand side of Eq. ([3]) . Backward scattering 
is synonymous with large momentum transfer, which in turn means corrections to RPA, which then translates to 
quadratic corrections in Eq.Q. Now we wish to compute the Green function of the Wigner crystal whose momentum 
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distribution has been computed in an earlier work^. There we considered electrons on a circle interacting via a long 
range interaction V(x) = —^\x\, where \x\ is the chord length. The ground state of this system was shown to be 
crystalline with lattice spacing l c = n/kp. It was shown that the momentum distribution of the Wigner crystal at 
zero temperature is given by, 

1 / "Wfl \ 1 



1 + e ^-* 2 \n F {p) + -[l + e ) (1 - n F (jp)) (13) 



where, w = This was derived using the general formula for the momentum distribution in terms of sea- 

bosons, 



( p ) + g-2 E q <^ +q/2 (q)A p + q/2 (q)>^ + (1 _ nF(p)) 1 ^ _ e -2 E q <A;_ q/2 (q)A p _ q/2 ( q) >^ (M) 



First we wish to generalize Eq. (|14p to finite temperature. This is important since we know that even for a nonin- 
teracting system, the momentum distribution at absolute zero is discontinuous at k — kp but is continuous at finite 
temperature. A naive approach that just performs thermodynamic averaging over the sea-boson occupation fails. We 
have to adopt a more subtle approach. One possibility is to express the sea-bosons in terms of the hydrodynamic 
variables and use the action in Eq.([5]). This can possibly be made to work out but is really not worth the effort. 
It is better to gain some intuition from this effort to guess the proper generalization. We propose the following 
generalization. 

n P = n F (p) ~ (l + A(p) e - 2 ^« A p + ^ (q) ^/ 2(q)>> ) + (l-n F (p)) i(l-A(p) e -^« A l-,M^-,M^ (15) 

Here A(p) contains the temperature information of the free theory only. The boson occupation in the exponent is 
defined as follows : <C ... ^>=< ... >p — < ... >p,o- That is, the difference between the interacting theory at finite 
temperature and free theory at finite temperature. To calculate the finite temperature sea-boson occupation we have 
to invoke the hydrodynamic description. We first express the sea-bosons in terms of hydrodynamic variables, namely 

the velocity potential and density. It is given as follows A k (q) = [A k (q), A k (q)] ("^[P-q — iXqj- Here we have to 

make sure that only the s-wave contributes, in other words we ignore complications caused by the dot product k.q 
and replace it by its extremum value. This correspondence automatically reproduces the RPA level identities also 
valid only in the s-wave sense : X)k(^k(q) + A^— = P-q an d J2k if]ff(^k(q) — A^— = ^q- We may also 
re-express these directly in terms of Fermi fields, p q = Z)k c ik+q/2 c k-q/2 and Xq = £ k sgn(k.q) c k - q /2 c k+q/2- 
This is beneficial since we may now easily compute the correlation functions of these operators at finite temperature. 
In general, < A k (q)A k (q) >= L4 k (q), A fc (q)] < p q/ 9_ q > + < A_ q A q > -j^j- For a Luttinger liquid the 

zero temperature theory is, 

H = E ^4(q)A k (q) + E Y[A{-q)A{q) + A\ q )A\-q)] = £ v\q\dt(q)d c (q) (16) 

k,q q q 

where, 



««> = , & r-^Y «i) ^ ah-,) 




where v = \] v 2 F - $ and [d c (q), 4(q)] = 1 and = A(q) + A^{-q) and X q = ^(Aiq) - A^(-q)) and 
Ml) = EkMQ) and [A(q),AHq)} = Isjk, p_ q = (d c (q) + 4(-q))c a and X q = (d c (q) - d\{-q))c b , c a = 

^((^"-(V) 1 ) and <% = vW((^) 4 + (^)")^- In other words < < >= 
\Ak(q) , A\(q)\(^ — l)]^fx- Substituting this into Eq. (TT3]) leads to a momentum distribution that has power law 
singularities with anomalous exponent 7 = ^ — 1. In general, we have to use the RPA- level action in Eq.© to 



6 



calculate the finite temperature expectation values namely, < p q p- q > and < JT q X_ q >. Let us now use this method 
to find the proper generalization of Eq. (fT3")) to finite temperature. First, we have to find A(p). By making contact 



with the free theory we find, A(p) = [sgn(k{ 
obtain, 



Tanh 



(P - e P ) 



When this procedure is implemented we 



A(p) e 



n F(p) + g ( 1 



Hp) e 



(1 - n F {p)) 



(19) 



Therefore the essential singularity remains even at finite temperature. Finally we wish to calculate the dynamical 
propagator of the Wigner crystal. It seems that in this case even Eq.© is not sufficient. We have to express the 
field operator in momentum space directly in terms of the sea-bosons. A partially correct formula was proposed in an 
earlier work&. Instead of searching for a rigorous approach, we directly make the following surmise for the propagators 
that is motivated by comparing with limiting cases. Define c p < = n F {p) c p and c Pi > = (1 — n F (p)) c p . Then, 



< c, 



(t )cp,>(t) >= (l-n F 



Cp)) \ ( 



1 - A(p) e 



-2£ q <^t (q,i)A p _ q/2 (q,i)>\ e -ie P (t-t) 



>E(\p\-k F ,t-t) (20) 



< c P ,<(*')ep,<W >= Mp) \ (l + A(p) e-^^l^M^^M^ 



-ia F {t-t ) 



E(\p\-k F ,t -t) 



(21) 



Here the envelope function has to be chosen with care. We make the following definition which we justify a posteriori 

r °° dwW(g,w)e-*"<*-*'> 



E(q,t-t) = 



$™duW{q,w) 



where W(q,tu) is the spectral weight, 



W(q, w) — Im 



(22) 



(23) 



e{q,u- i5) i 

In one dimension, the spectral weight is a delta function at the collective mode. Hence in one dimension, 

E{q,t-t) = e -^(D(t-t') (24) 

In more than one dimension we have both the particle-hole mode and collective mode. For example in the case of 
the jellium, the collective mode contribution occurs at the plasma frequency and this leads to a rapidly oscillating 
contribution which may be ignored. The important contribution comes from the particle-hole mode, which even 
though is not infinitely long-lived, makes a significant contribution. For the noninteracting theory in any number of 
dimensions this prescription gives us, 



E free (q,t-t)=e- i °*M*-t) 



(25) 



One may then use the Kubo-Martin-Schwinger (KMS) boundary conditions to evaluate < c p (t)c^ p (t ) > from these 
propagators. First we observe that this prescription reproduces the dynamical Green function of the free theory 
in any number of dimensions both at finite temperature and at zero temperature. We may use these to compute 
the dynamical propagator of the Luttinger liquid and see if these results agree with those of the more traditional 
approaches. Using the traditional method, 



< i/j^x ,t')ip R (0,t) >= e 



— „—ieF(t—t ) _ 



2ni [x - v(t' - t)] 
Taking the Fourier transform with respect to x we obtain, 



1 



= p -ie F {t-t) -i{\p\-kp)v{t-t) 1 



1 + sgn(k F - \p\) 



A\x - v(t' - t)\ 

fc F -bir T 



A 



(26) 



(27) 



This is completely identical to the result using Eq. (|2"D|) and Eq. (l2"l"T) since uj(q) = v\q\. This approach may seem quite 
ad-hoc but it would be very surprising indeed if a theory that reproduces all dynamical aspects of the free theory 
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both at finite temperature and zero temperature correctly in any number of dimensions and also able to reproduce 
the very nontrivial dynamical propagator of the Luttinger liquid correctly in one dimension is not valid in general. 
Besides, for the jellium, the clever choice of the envelope function ensures that one-particle Green function is gapless 
even though the collective mode is gapped. Assuming that this is valid in general we may write down the dynamical 
propagator of the Wigner crystal as follows. 

< cl(t')c p (t) >= I (l + A(p) e " W"^^ n F (p)e-"-( t - t ')e^( t -*') 

+ ~ (l - A(p) e -^h coth ^\ (i _ Uir ( F )) e -^(*-*') e -^o(*-*') (28) 



< c 



>= 



1 - A(p) e 



oth( 



n F (p)e- ieF(t - t ) e iu °^- t ^ 



+ + A(p) e ^ c ° t,l( (l - npip^e-^-^e-^-^ (29) 

Note that the above formulas are independent of arbitrarily chosen momentum cutoffs and depend only on the 
microscopic parameters present in the original hamiltonian with parabolic dispersion. This in contrast with the 
conventional approaches to bosonizationi where such cutoffs are mandated by the formalism. In this work, Schulz 
remarks that the time dependent and temperature dependent formulas for the propagators are complicated in his 
formalism and hence not too illuminating. This is in contrast with our 'momentum space bosonization' approach 
where the above general formulas are not only simple but also illuminating. 



IV. CONCLUSIONS 



To conclude, we have computed the dynamical Green function of the Wigner crystal in one dimension whose 
momentum distribution exhibits essential singularities in momentum space. We have generalized the momentum 
distribution to finite temperature. We have written down a general formula for the field operator without any Klein 
factors in one dimension in terms of currents and densities that does not involve momentum cutoffs and applies directly 
to the Fermi gas with a parabolic dispersion rather than to its caricature namely the Luttinger model. Lastly we 
have summarized all the developments in the subject made by the author and his collaborators in order to facilitate 
further developments. 
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